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Abstract 

The deformation barriers standing in the quasi-molecular shape 
path have been determined in the actinide region within a macroscopic- 
microscopic energy derived from a generalized liquid drop model, the 
algebraic droplet model shell corrections and analytic expressions for 
the pairing energies. Double and triple-humped fission barriers ap- 
■ pear. The second barrier corresponds to the transition from one-body 

shapes to two touching ellipsoids. The third minimum and third peak, 
when they exist, come from shell rearrangements in the deformed frag- 
ment. The shape of the other almost magic one is close to the sphere. 
The barrier heights agree with the experimental results, the energy of 
the second minimum being a little too high. The predicted half-lives 
follow the experimental data trend. 
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O ■ 1 Introduction 
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The possibility of transmutation of nuclear waste and of production of en- 
ergy by accelerator-driven systems is under consideration. The knowledge 
of all the nuclear reactions which constitute a non negligible part of the re- 
action cross section is needed pQ. Different codes are under construction or 
improvement (Fluka, Gnash, Talys,..) and accurate potential barriers must 
be calculated rapidly, particularly in the actinide region, to predict or firstly 
to reproduce the mass and charge distributions governing the fission cross 
sections [2|. Furthermore, new precise measurements renew also interest in 
investigating the multiple-humped barriers of the actinide nuclei and heaviest 
elements. The analysis of the fission probability and of the angular distribu- 
tion of the fission fragments support the presence of hyperdeformed states 
in a deep third well in several Th and U isotopes |S1 lU E] confirming the 
pioneering work of Blons et al jH] in 231 > 233 Th. The observed strongly en- 
hanced low energy a decay in some heavy actinide nuclei is also explained 
by transition from a third hyperdeformed minimum and the possibility that 
the third minimum is the true ground state of very heavy and perhaps super- 
heavy nuclei is even also advocated [7j. In medium mass nuclei some signs of 
hyperdeformed rotational bands have been found, but no discrete HD level 
has been identified [HJ E] • 
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By adding at the macroscopic liquid drop model energy of elongated one- 
body shapes an oscillatory microscopic contribution, the Strutinsky's method 
PE] generated double-humped barriers allowing to predict and explain the 
fission isomer characteristics. Myers and Swiatecki P™1J proposed analytic 
formulae in the same objective of taking into account the shell and pairing 
energies. Later on, the asymmetric two-center shell model J2] , Hartree-Fock- 
Bogoliubov :Wl and relativistic mean field theories jHj have also allowed to 
reproduce double-humped barriers. The heights of the inner and asymmetric 
outer fission barriers are almost constant (5-6 MeV) from Th to Am isotopes 
[TH| llfij. It is a severe test for the theoretical models. In the actinide region 
the third hyperdeformed minimum was predicted by Moller et al and by 
recent theoretical approaches [T8J EH I2D1 • 

It has been previously shown within a Generalized Liquid Drop Model 
taking into account both the proximity energy between close opposite sur- 
faces, the asymmetry and an accurate radius that most of the symmetric 
and asymmetric fission [2H 123 , <x and light nucleus emissions [23 121] and 
super and highly deformed state data j2S] can also be reproduced in the 
fusiolike shape path. The purpose of this work is, within this GLDM, to 
go beyond the two separated sphere approximation by taking into account 
the ellipsoidal deformations of the two different fission fragments and their 
associated shell and pairing energies, investigating all the possible mass and 
charge asymmetries. 



2 Potential energy 

The total energy of a deformed nucleus is the sum of the macroscopic GLDM 
energy, the shell correction energy and the pairing energy. The GLDM energy 
is given by [2E] 

E = E v + E s + E c + E prox + E rot , (1) 

where the different terms are respectively the volume, surface, Coulomb, 
nuclear proximity and rotational energies. 

All along the deformation path the nuclear proximity energy term E prox 
allows to take into account the effects of the attractive nuclear forces between 
nucleons facing each other across a neck in the case of a deformed one-body 
shape or across a gap in the case of two separated fragments. This is not 
a small correction in the quasi-molecular shape path. For example, at the 
contact point between two spherical Kr and Ba nuclei the proximity energy 
reaches —43 MeV. 

E wox {r) = 2 7 1 $ [D{r, h)/b] 2nhdh. (2) 
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r is the distance between the mass centres, h is the transverse distance 
varying from the neck radius to the height of the neck border. D is the 
distance between the opposite surfaces and b the surface width. <3> is the 
proximity function and 7 the surface parameter. 

2.1 One-body-shapes 

For one-body shapes, the first three contributions are given by 

E v = -15.494(1 - 1.8I 2 )A MeV, (3) 

E s = 17.9439(1 - 2.6/ 2 U 2/3 -4^ MeV, (4) 

AttRq 

E c = O.Qe 2 (Z 2 /R )B c . (5) 

B c is the Coulomb shape dependent function, S is the surface and / is the 
relative neutron excess 

B c = 0.5 J (V(9)/V Q ) (R(6)/R ) 3 sin 6d6, (6) 

where V(6) is the electrostatic potential at the surface and Vo the surface 
potential of the sphere. 

The radius R of the compound nucleus is defined as: 

R = (1.28A 1/3 - 0.76 + 0.8A~ 1/3 ) fm. (7) 

which leads, for example, to Rq = 5.3 fm and r = 1.15 fm for 98 Zr and 
Rq = 7.5 fm and r = 1.18 fm for 255 Fm. The radius of the two fragments is 
calculated assuming volume conservation. 

As in previous works, the one-body shape sequence is described within 
two joined elliptic lemniscatoids which allow to simulate the development of 
a deep neck in compact and little elongated shapes with spherical ends (see 
Fig. HJ. The proximity energy is maximized in this deformation path. 

2.2 Two separated ellipsoids 

For two-body shapes, the coaxial ellipsoidal deformations have been consid- 
ered j2Z] (see Fig. |2J). The system configuration depends on two parameters 
: the ratios (i = 1, 2) between the transverse semi-axis aj and the radial 
semi-axis q of the two different fragments. 

di = Ris}^ 3 and Ci = RiS { 2 ^ 3 . (8) 
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Figure 1: Selected shape sequence to simulate the one-body shape evolution. 




Figure 2: Two coaxial ellipsoid configuration describing the two-body shape 
part of the fission barrier. The fission axis is the common axis of revolution. 



The prolate deformation is characterized by s < 1 and the related eccen- 
tricity is written as e 2 = 1 — s 2 while in the oblate case s > 1 and e 2 = 1 — s~ 2 . 
The volume and surface energies are Ey 12 = Ey 1 + Ey 2 and Es 12 = Es 1 + E$ 2 . 
In the prolate case, the relative surface energy reads 




.2N1/3 r 



sin 1 (e i ) 



2 



1 + 



e,(l-e 2 )V2 



(9) 



and in the oblate case 



(1 + 6 2 )V3 



1 + 



ln(e t + (l + e 2 y/ 2 ) 
e t (l + e 2 y/ 2 




- 1. 



(10) 
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The Coulomb self-energy of the spheroid i is 



Se 2 Z 2 B a 



(11) 
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The relative self-energy is, in the prolate case 



B Ci = [ - 6J In—^- (12) 
° 2e, l-e, : K J 



and, in the oblate case 



i? C4 = 1 - J tan- 1 ^. (13) 
The Coulomb interaction energy between the two fragments reads 

E c ,int = [s(Xi) + s(X 2 ) - 1 + S(A 1; A 2 )] A, 2 = (14) 



r being, as before, the distance between the two mass centres. 
In the prolate case, s(Aj) is expressed as 



= irr (15) 



while, for the oblate shapes 

3 1 1 

2 Ui uf 2u. 



s (Xi) = -(— + —)tan 1 u i - — a; 2 = -A 2 . (16) 



^(Ai, A2) can be represented in the form of a two- fold summation 
5(Al ' A2) = § £ (2j + l) 3 (2j + 3) (2fc + l) 3 (2fc + 3) %Wk)\ A W ' (17) 

3 Shell energy 

The shell corrections for a deformed nucleus have been determined within 
the algebraic formulae given in the Droplet Model [2H] with slightly different 
values of the parameters. 

E sheU = E:tjra-^o 2 )e-°- 9e2 . m 

The shell corrections for a spherical nucleus are 

E s s p h h e a e = 5.8 [(F(N) + F(Z))/(0.5A) 2 / 3 - 0.26A 1 / 3 ] MeV, (19) 
where, for M»_i < X < M t) 

F{X) = qi (X - M^) - 0.6(X 5 / 3 - M?H). (20) 
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Mi are the magic numbers and 



qi = 0.6(Mr - M^i)/(Mi - Mi-t). (21) 

The selected highest proton magic number is 114 while, for the two highest 
neutron magic numbers, the values 126 and 184 have been retained. 

6 2 = (5R) 2 /a 2 . (22) 

The distortion 9a is the root mean square of the deviation of the nuclear 
surface from the sphere, a quantity which incorporates indiscriminately all 
types of deformation. The range a has been chosen to be 0.32ro- 
For the two-body shapes, the total shell energy is the sum of the shell cor- 
rections for each deformed fragment. 

4 Pairing energy 

The pairing energy has been calculated with the following expressions. 
For odd Z, odd N and N=Z nuclei 

E Pairing = 4.8/iV 1 / 3 + 4.8/Z 1/3 - 6.6/ A 2/3 + 30/ A. (23) 

For odd Z, odd N and N 7^ Z nuclei 

Ep ai rin g = 4.8/iV 1 / 3 + 4.8/Z 1 / 3 - 6.6/ A 2 / 3 . (24) 

For odd Z, even N nuclei 

Epairing = 4.8/Z 1 / 3 . (25) 

For even Z, odd N nuclei 

Ep ai ring = 4.8/^/3. ( 2 6) 

For even Z, even N nuclei 

Epairing 0. (^7) 

5 Potential barriers 

The dependence of the potential barriers on the shape sequence and introduc- 
tion of the microscopic corrections is displayed in Fig. El for an asymmetric 
fission path of the 230 Th nucleus. The shell effects generate the slightly de- 
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Figure 3: Asymmetric fission barrier of a 230 Th nucleus emitting a doubly 
magic nucleus 132 Sn. The dashed and dashed-dotted curves give the energy 
within the two-sphere approximation for the two-body shapes without and 
with shell corrections around the original sphere while the dotted and solid 
lines include the ellipsoidal deformations without and with shell energies, r 
is the distance between mass centres. 
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formed ground state and contribute to the formation of the first peak. The 
proximity energy flattens the potential energy curve and will explain the 
formation of a deep second minimum for heaviest nuclei. The shell effects 
are attenuated at these large deformations. In the exit channel correspond- 
ing to the two-sphere approximation the top of the barrier (r = 12.6 fm on 
this example) is reached after the rupture of the matter bridge between the 
two spherical fragments (r = 11.4 fm). Then, the top corresponds to two 
separated spherical fragments maintained in unstable equilibrium by the bal- 
ance between the attractive nuclear forces and the repulsive Coulomb ones. 
In this path, the introduction of the shell and pairing effects for two-body 
shapes is not sufficient to reproduce accurately the experimental data on the 
fission barrier heights of actinide nuclei. The transition between one-body 
and two-body shapes is less smooth when the ellipsoidal deformations of the 
fragments and the proximity energy are taken into account. It corresponds 
to the passage (at r = 11 fm for 230 Th) from a one-body shape with spherical 
ends and a deep neck to two touching ellipsoidal fragments, one or both of 
them being slightly oblate. The barrier height is reduced by several MeV. 
The introduction of the shell effects still lowers the second peak and shifts it 
to an inner position (r = 10.3 fm here). It even leads to a third minimum 
and third peak in this asymmetric decay path. A plateau appears also at 
larger distances around 10 MeV below the ground state. It is due to the end 
of the contact between the two fragments and the persistence of the prolate 
deformation of the lightest fragment. The end of the plateau corresponds 
to the rapid transition from prolate to oblate shapes for the non-magical 
fragment and the vanishing of the proximity energy. This second fragment 
returns to a prolate shape when the interaction Coulomb energy is smaller. 
The potential barriers for the 232,235,238^ 238,240,243p U; 243,244 Am> 243 Cni) 

250 Bk and 250 Cf nuclei are shown in Fig. @] to |7| It is important to men- 
tion that to obtain these barriers the only input parameters are Ai,Zi,A 2 
and Z 2 and that the calculations are very rapid on an usual computer and, 
consequently, can be integrated in a more general and complex code. 

For a given mass asymmetry, the charge asymmetry which minimizes the 
deformation energy has been selected. The proximity energy and the at- 
tenuated microscopic effects are responsible for the formation of a second 
one-body shape minimum. The heights of the two peaks generally increase 
with the asymmetry but the shell and pairing corrections induce strong vari- 
ations from this global behaviour. Their main effect is to favour, for the U, 
Pu, Am and Cm isotopes, an asymmetric path where one fragment is close to 
the doubly magic number 5Q 2 Sn nucleus, and, consequently, keeps an almost 
spherical shape. This effect is less pronounced for 243 Cm and 250 Cf since for 
nuclei with Z ~ 100 the symmetric fission gives fragments with a charge of 
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around 50. For these nuclei a symmetric and an asymmetric exit channels 
are compared on the lower part of the Figures El and [7| A third minimum and 
third peak appear in the asymmetric decay path. There is no third barrier in 
the symmetric deformation path. The origin of the third peak is investigated 
in section VI. 

The calculated and experimental energies of the maxima and minima of 
the fission barriers are compared in table 1. The choice of the most probable 
fission path is difficult for some isotopes since there is a true degenerescence 
in energy between several paths of the mult i- dimensional potential surfaces, 
particularly for the heaviest elements where the symmetric path seems more 
probable. There is a very good agreement between the experimental and 
theoretical heights E a and E b of the two peaks. The predicted value of the 
second minimum energy is a little too high. Two reasons may be advanced 
to explain this difference. Firstly, the shell energy may be underestimated 
at this large deformation. Secondly, the shape sequence may introduce too 
rapidly the asymmetry, which would lead to an undervaluation of the proxim- 
ity energy. The still sparse data for the third barrier are correctly reproduced. 
For the heaviest nuclei the external barrier disappears since the attractive 
proximity forces can no more compensate for the repulsive Coulomb forces. 



9 




r (fm) A (heaviest) 




6 7 8 9 10 11 12 120 125 130 135 140 145 150 155 

r (fm) A (heaviest) 



Figure 4: On the left, multiple- humped fission barriers in the mentioned 
asymmetric fission path for 232 > 235 > 238 u. Q n ^he right, inner (full circles) and 
outer (crosses) fission barrier heights as a function of the mass of the heaviest 
fragment. 
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Figure 5: On the left, multiple-humped fission barriers in the mentioned 
asymmetric fission path for 2 38,24o,243p u ^ Q n ^ e right, inner (full circles) and 
outer (crosses) fission barrier heights as a function of the mass of the heaviest 
fragment. 
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Figure 6: Fission barriers for 243,244 Am and 243 Cm. The inner and outer 
fission barrier heights are given, on the right, respectively by the full circles 
and crosses. 
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Figure 7: Fission barriers for 250 Bk and 250 Cf. The inner and outer fission 
barrier heights are given, on the right, respectively by the full circles and 
crosses. 
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Table 1. Comparison between theoretical and experimental [T51 1^ \W\ 14"] 

barrier characteristics for actinide nuclei. E a , Ej, and E c are the first, second 
and third peak heights while En and Em are the energies of the second and 
third potential minima relatively to the ground state energy (in MeV). 

6 Third barrier 

The origin of the existence of the third well in the asymmetric decay path is 
examined now (see Fig. |HJ). The dashed line represents the potential for two 
touching ellipsoids when the one-body shape is still energetically favoured. 
The second peak (but first on the figure) corresponds to the point where these 
touching ellipsoids begin to give the lowest energy. The heaviest fragment 
is a magic nucleus. It therefore preserves its almost spherical shape. The 
non magic fragment was born in an oblate shape (s ~ 1.4), due to the small 
distance between the mass centres at this step. When this distance increases, 
the ratio s decreases, because of the proximity energy which tends to keep 
close the two tips of the fragments. Thus, the lightest fragment remaining 
in contact with the other spherical fragment approaches the spherical shape 
and its shell energy increases to reach a maximum which is at the origin 
of the third peak and which corresponds to two touching different spheres. 
Before reaching this third peak a third minimum appears. Its shape is hy- 
perdeformed and asymmetric in agreement with the experimental data jlj. 
Later on, the proximity forces maintain the two fragments in contact and 
the shape of the smallest one evolves to prolate shapes (s < 1) and the shell 
corrections decrease. The third barrier appears only in the asymmetric de- 
cay path and for some specific nuclei. In the symmetric mass exit path, the 
proximity and Coulomb energies counterbalance the smallest shell effects and 
induce an asymmetric shape, the two fragments remain in contact but one 
fragment is oblate while the other one is prolate. With increasing distance 
between the mass centres the two nuclei become prolate. 

The dependence of the fission barrier heights and profiles on the asym- 
metry for the 231 ' 233 Th and 234 ' 236 U nuclei, for which experimental data on 
the third barrier exist, are given in Fig. El and El The position of the second 
peak in the symmetric decay path corresponds to the position of the third 
peak in the asymmetric deformation path. Clearly the magicity of some Sn 
isotopes plays the main role. 
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Figure 8: Fission barriers, shell energies and ratio of the semi-axes of the two 
ellipsoidal fragments for an asymmetric decay channel and the symmetric 
one for 236 U. On the lowest part, the fission barrier is given by the solid line. 



7 Half- lives 

Within this asymmetric fission model the decay constant is simply given by 
A = vqP. The assault frequency u has been taken as z/ = 10 20 The 
barrier penetrability P is calculated within the action integral 

2 r° ut i 

P = exp [ J2B{r){E{r)-E g . s )dr}. (28) 

Tl Jri„ 



The inertia B(r) is related to the reduced mass by 

B(r) = n[l + 24exp(-3.25(r - R sph )/R )} (29) 

where R sp h is the distance between the mass centres of the future fragments 
in the initial sphere, R sp h/Ro = 0.75 in the symmetric case. 
The inertia depends on the internal structure of the system. A large amount 
of internal reorganization occuring at level crossings raises the inertia. For 
shapes near the ground state the inertia is expected to be considerably above 
the irrotational flow value. For shapes remaining a long time highly deformed 
the reduced mass is reached asymptotically. The selected phenomenological 
inertia B(r) adjusted to reproduce the experimental data is compared with 
other previous semi-empirical inertia in Fig. ^2 Its value is slightly higher 
at the beginning of the fission process. It reaches more rapidly the reduced 
mass value for elongated shapes. 
The partial half-life is finally obtained by T±/ 2 = 
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Figure 9: Fission barriers for 231 > 233 Th and two different decay channels. The 
inner and outer fission barrier heights are given, on the right, respectively by 
the full circles and crosses. 
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Figure 10: Fission barriers for 234 > 236 U. The inner and outer fission barrier 
heights are given, on the right, respectively by the full circles and crosses. 
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Figure 11: Comparison between different selected semi-empirical inertia B(r). 
The dashed line corresponds to the inertia proposed in ref. |2H| while the 
dashed and dotted curve and dotted line give respectively the inertia pro- 
posed in ref. [3U] for elongated shapes and for compact and creviced shapes. 
The values obtained from the formula (29) are given by the solid line. 
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The experimental spontaneous fission half-lives and theoretical predic- 
tions for the supposed most probable exit channels are compared in Table 
2. The half-lives corresponding to several paths are given for 243 Cm, 250 Cf, 
256p m anc j 256]\j a Except for the ligthest U isotopes, there is a correct agree- 
ment with the experimental data on 20 orders of magnitude. The half-lives 
vary regularly for close exit channels of a same nucleus. 



Reaction 




Ti/2,theo(s) 


232/7 134m ,98 y r 
92 u 52 c '40 


2.5 x 10 21 


3.6 x 10 16 


234^ IfSn+lfMo 


4.7 x 10 23 


8 x 10 19 


IfU^ IfSn+l^Mo 


3.1 x 10 26 


7.7 x 10 23 


IfU^ l 5 fSn+l^Mo 


7.8 x 10 23 


1.0 x 10 22 


238 f/-> IfSn+lfMo 


2.6 x 10 23 


5.3 x 10 22 


IfPu -> $°Sn +\f Ru 


1.5 x 10 18 


2.6 x 10 19 


IfPu -> f ° Sn +™ Ru 


2.5 x 10 23 


4.8 x 10 22 


IfPu ™>Sn +™ Ru 


3.7 x 10 18 


4.8 x 10 19 


™*Am^ ^fSb+^Ru 


6.3 x 10 21 


1.1 x 10 23 


mCm^ ™Sn+™Pd 


1.7 x 10 19 


3 x 10 21 


Q^Cl7l ► 4g 2 Cc? ~t"48^ Cd 


1.7 x 10 19 


1.6 x 10 18 


If Cm -> l^Sn Pd 


4.4 x 10 19 


3 x 10 20 


IfCm^ IfSn+^Pd 


1.3 x 10 14 


7.7 x 10 15 


250^ f __. 125 r„ ,125 r„ 
98 °J ^ 49 J n +49 J n 


5.2 x 10 11 


1.9 x 10 9 


250^/ , 132m ,118 pj 
98 °J ^ 52 1 e +46 ra 


5.2 x 10 11 


1.2 x 10 10 


250^ f . 140^< „ ,110 rp 
98 ° / ~^ 55 ° S +43 J C 


5.2 x 10 11 


4.9 x 10 11 


2 f Es -> g»Sn +if In 


8.4 x 10 10 


8 x 10 9 




1.0 x 10 4 


45 


ggFm - ^ + 135 I 


1.0 x 10 4 


82 


f&No - g» £6 + 5 f 56 


110 


0.9 x 10~ 2 


ggJVo- %*Pd+%>Ba 


110 


0.3 x IO" 1 



Table 2. Comparison between experimental PB] and theoretical 
spontaneous fission half-lives of actinide nuclei. 

8 Summary and conclusion 

Double and triple-humped fission barriers appear for the actinide elements 
in the quasi-molecular shape path within a macroscopic-microscopic defor- 
mation energy derived from a generalized liquid drop model and analytic 
expressions for the shell and pairing energies. The second peak corresponds 
to the transition from one-body shapes to two touching ellipsoids. The third 
barrier appears only in the asymmetric decay path and for some specific nu- 
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clei. Then, the heaviest fragment is almost a magic nucleus and it preserves 
its shape close to the sphere. The other fragment evolves from an oblate 
ellipsoid to a prolate one and the third peak corresponds to the maximum of 
the shell effects in the non magic fragment and, consequently, to two touching 
different spheres. The barrier heights agree with the experimental results for 
the double and triple-humped fission barriers, the energy of the second mini- 
mum being a little too high. The predicted half-lives follow the experimental 
data trend. 
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